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One-band Hubbard model with hopping parameter t and Coulomb repulsion U is considered at 
half filling. By means of the Schwinger bosons and slave Fermions representation of the electron 
operators and integrating out the spin-singlet Fermi fields an effective Heisenberg model with an- 
tiferromagnetic exchange constant is obtained for vectors which identifies the local orientation of 
the spin of the itinerant electrons. The amplitude of the spin vectors is an effective spin of the 
itinerant electrons accounting for the fact that some sites, in the ground state, are doubly occu- 
pied or empty. Accounting adequately for the magnon-magnon interaction the Neel temperature is 
calculated. When the ratio jj is small enough (77 < 0.09) the effective model describes a system 
of localized electrons. Increasing the ratio increases the density of doubly occupied states which in 
turn decreases the effective spin and Neel temperature. The phase diagram in plane of temperature 
^ff- and parameter jj is presented. The quantum critical point (Tn = 0) is reached at jj = 0.9. The 
magnons in the paramagnetic phase are studied and the contribution of the magnons' fluctuations 
to the heat capacity is calculated. At Neel temperature the heat capacity has a peak which is 
suppressed when the system approaches quantum critical point. 
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PACS numbers: 64.70.Tg, 71.10.Fd, 75.50.Ec, 75.40.Cx 



I. INTRODUCTION 

Quantum phase transitions arise in many-body sys- 
tems because of competing interactions that support 
different ground states. At quantum critical point 
(QCP)thc matter undergoes a transition from one phase 
to another at zero temperature. A nonthermal control 
parameter, such as pressure, drives the system to QCP. 
The typical temperature-pressure phase diagrams ob- 
served in the heavy- fermion materials CePd 2 Si 2 , Celn^, 
CeRh 2 Si 2 , CeCu 2 Si 2 and CeRhIn 5 show that at 

ambient pressure the compounds order into antifcrromag- 
nets below the Neel temperature TV. Applying pressure 
reduces TV monotonically. The QCP is the critical pres- 
sure at which the Neel temperature TV = 0. 

The quantum critical behavior has been extensively 
studied for many years. Many books @,||| , review articles 
(9I-HH and papers investigate systematically the magnetic 
quantum critical point. 

The magnetism of cerium based compounds is deter- 
mined by the 4/ electrons of Ce 3+ ions. The strong spin- 
orbit coupling splits the 4/ electrons into j ' = | and 
] = \ multiplcts , where j is the total angular momen- 
tum. Only the sextuplet effectively contributes to the low 
energy excitations. It is further split into T7 doublet and 
Ts quadruplet due to crystal electric field. For isotropic 
systems like Celnj,, the energy level of T7 is lowest. The 

eigenstates are \Tj± >= ± | > — y§| T § >, where 
"+" and "-" denote up and down "pseudo-spins" respec- 
tively. The three-dimensional (3D) one-band Hubbard 
model is the simplest model of itinerant magnetism of 
the isotropic cerium based compounds. Although the hy- 
bridization with Inbp electronic states my be important, 
here it is considered as a renormalization of the hopping 
amplitude. 
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The ground state of the 3D Hubbard model on a sim- 
ple cubic lattice, at half filling has antiferromagnetic long 
range order for all positive values of the onsite Coulomb 
repulsion. There have been many attempts to calculate 
Neel temperature TV using quantum Monte Carlo(QMC) 
simulations Il^|-{l9| , variational methods [20. 
Fock theory [22| , strong coupling expansions 
namical mean field theory (DMFT) [25l - [27j . Calculations 
beyond the dynamical mean field theory: the dynamical 
cluster approximation (DCA) [2a |, cluster generalization 
of the DMFT [29| and dynamical vertex approximation 
(DTA) [3(1 [HJ have been proposed. All these studies do 
not show any trace of quantum critical point. 

Here we focus attention on the quantum critical be- 
havior in itinerant antifcrromagncts. The increasing of 
the double occupancy in 31? one-band Hubbard model 
at half filling pushes the system to the quantum criti- 
cality. Increasing the ratio jj, where t is the hopping 
parameter and U is the Coulomb repulsion, increases the 
density of doubly occupied states which in turn decreases 
the effective spin and Neel temperature. The Quantum 
Critical Point is a state with a critically high density of 
the doubly occupied state. 

The physics of the antiferromagnets near the Neel tem- 
perature is dominated by the magnons' fluctuations and 
it is important to account for them in the best way. We 
employ a technique of calculation, which captures the es- 
sentials of the magnons' fluctuations in the theory, and 
for 2D systems one obtains zero Neel temperature, in 
accordance with Mermin- Wagner theorem |33|. 

The paper is organized as follows. In Sec. II, starting 
from one-band Hubbard model at half filling ,with hop- 
ping parameter t and Coulomb repulsion U, we derive an 
effective Heisenberg-like model, with antiferromagnetic 
exchange constant, in terms of the vector describing the 
local orientations of the magnetization. The transversal 
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fluctuations of the vector are the magnons in the the- 
ory. The amplitude m of the spin vectors is an effective 
spin of the itinerant electrons accounting for the fact that 
some sites, in the ground state, are doubly occupied or 
empty. This is a base for Neel temperature calculation. 
Section III is devoted to T/v/J — m and Tn/U — t/U 
phase diagrams of the model. The paramagnetic phase 
of itinerant antifcrromagnets is explored in Section IV. 
The calculations of the specific heat for different values 
of the control parameter t/U are presented in the Section 
V. A summary in Sec. VI concludes the paper. 

II. EFFECTIVE MODEL 

We consider a theory with Hamiltonian 

/i = -i E { c t* c j<y + h - c ) + U Y1 n ^ n H -mE" 1 w 
(ij) 1 1 

where cf a and Ci a (a =f, ]/) are creation and annihilation 
operators for spin-1/2 Fermi operators of itinerant elec- 
trons, riia- = cl a Ci„, m = + ni, t > is the hopping 
parameter, U > is the the Coulomb repulsion and /i 
is the chemical potential. The sums are over all sites of 
a three-dimensional cubic lattice, and denotes the 
sum over the nearest neighbors. 

We represent the Fermi operators, the spin of the itin- 
erant electrons 

•< = ^E c -^ c — ( 2 ) 

aa' 

where (t x , t v ,t z ) are Pauli matrices, and the density op- 
erators riic in terms of the Schwingcr bosons (<£j )Cr , ft a) 
and slave Fcrmions (h i} ht 7 d i} dt). The Bose fields are 
doublets (a = 1,2) without charge, while Fermions are 
spinless with charges 1 (di) and -1 (hi). 

Ci-J- = htfn + tpt 2 dii c 4 = h tVi2 ~ ftidh 

ni = 1 - hfhi +dt~di, s\ = - ^ fta T aa' <Pio' , 

aa' 

c t c n c t c a = d t d i (3) 

ftifii + <Pi2<Pi2 + dtdi + hfhi = 1 (4) 

To solve the constraint (Eqf?]), one makes a c hang e of 
variables, introducing Bose doublets £ ?(T and Cit |34| 

Cia = <Pia (1 - hfhi - dtdi) 2 , 

C = (5) 

where the new fields satisfy the constraint C 8 >Cia = 1- In 
terms of the new fields the spin vectors of the itinerant 
electrons Eq. (J2j) have the form 

"i= \ E C 1 - h tfH - 4di\ (6) 



When, in the ground state, the lattice site is empty, the 
operator identity hthi = 1 is true. When the lattice 
site is doubly occupied, dtdi = 1. Hence, when the lat- 
tice site is empty or doubly occupied the spin on this 
site is zero. When the lattice site is neither empty nor 
doubly occupied (hthi = dtdi = 0), the spin equals 
Sj = l/2rij, where the unit vector 

<-ECC(»' (nf = l) (7) 

aa 1 

identifies the local orientation of the spin of the itinerant 
electron. 

The Hamiltonian Eq.([T]), rewritten in terms of Bose 
fields Eq.© and slave Fermions, adopts the form 

h = -tJ2Mdi-hthi)(lQ° 

+ (d+ht - 4hf) (CnOa ~ CiaOO + h - c ] 

x (l - hthi - dtdi) * (1 - h+hj - d+dj) * 

+ U^4di -^(1- hfhi + dttk), (8) 

i i 

An important advantage of working with Schwinger 
bosons and slave Fermions is the fact that Hubbard term 
is in a diagonal form. The fcrmion-fcrmion and fcrmion- 
boson interactions are included in the hopping term. One 
treats them as a perturbation. To proceed we approxi- 
mate the hopping term of the Hamiltonian Eq.® setting 

(l — hthi — dt di) 2 ~ 1 and keep only the quadratic, 
with respect to Fermions, terms. This means that the 
averaging in the subspace of the Fermions is performed 
in one fermion-loop approximation. Further, we repre- 
sent the resulting Hamiltonian as a sum of two terms 

h = ho + h int , (9) 

where 

ho= - * E ( d t di - h i hi + h - c ) + u E d t ck 

{ij) " 

- ^(l-htlH + dtdi), (10) 

i 

is the Hamiltonian of the free d and h fermions, and 
hint = - * £ [(dfdi - hfhi) (CfQ* - 1) (11) 

{«> 

+ ( d t h t - d t h t) (CflC;a ~ CeOi) + h.c] 

is the Hamiltonian of boson-fermion interaction. 

The ground state of the system, without accounting for 
the spin fluctuations, is determined by the free-fermion 
Hamiltonian ho and is labeled by the density of electrons 

n = 1- < hfhi > + < dtdi > (12) 
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(see equation ([3]) ) and the " effective spin" of the electron 
m = i (1- < h+hi >-< d+d t >) . (13) 
At half filling 



< hfhi >=< dfck > 



(14) 



To solve this equation, for all values of the parameters 
U and t, one sets the chemical potential \i = U/2. Uti- 
lizing this representation of \i we calculate the effective 
spin "m" as a function of the ration t/U. The result is 
depicted in figure (JTJ. 

Let us introduce the vector, 



1 Cia T aa'Ci<y' 



(15) 



Then, the spin- vector of itinerant electrons Eq. (|6j can be 
written in the form 



1 

2m 



M, (1 - h+ hi - d+ di) 



(16) 



where the vector Mj identifies the local orientation of 
the spin of the itinerant electrons. The contribution of 
itinerant electrons to the total magnetization is < sf >. 
Accounting for the definition of m (see Eq fi3|) . one ob- 
tains < sf >=< Mf >. 

The Hamiltonian is quadratic with respect to the 
Fermions di , df and hi , hf , and one can average in the 
subspace of these Fermions (to integrate them out in the 
path integral approach). As a result, one obtains an ef- 
fective model for vectors M,, which identifies the local 
orientation of the spin of the itinerant electrons, with 
Hamiltonian 



(17) 



The effective exchange constant J is calculated in the one 
loop approximation and in the limit when the frequency 
and the wave vector are small. At zero temperature, one 
obtains 



J = 

t 1 

2t 2 1 

3m 2 UN 



(18) 



E C0S M [°(-4) + o(-4)] 



k V=l 
/ 3 



E E^H [l-e(-e h k )-9(-ei)} 



and the wave vector k runs over the first Brillouin zone 
of a cubic lattice. 

The exchange constant J and the effective spin m are 
function of the ratio t/U. At half filling the exchange 
constant J is positive and the model (fT7)) is an effective 
model of itinerant antiferromagnetism. The functions 
m(t/U) and J(t/U)/U are depicted in figure ([J). At half 
filling the density of doubly occupied states < d + d > is 
equal to the density of empty states < h + h >. Increasing 
the ratio t/U increases the density of doubly occupied 
states which in turn decreases the effective spin of the 
system (see equation ([13])). 



m 



♦ A 



J/U 



3,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,C 

t/U 



FIG. 1: (Color online) The effective spin of the system m 
as a function of the ratio t/U-black triangles (left scale). Di- 
mensionless exchange constant J/U as a function of t/U-red 
rhombuses (right scale). 



III. PHASE DIAGRAM 

We are going to study the antiferromagnetic phase of 
the model Eq. (fl"7)) with J > 0. To proceed, one uses 
the Holstein-Primakoff representation of the spin vectors 



Mj(ai~, dj), where aj 



dj are Bosc fields. 



Mf 



Mj 



M] 



cos 2 y y 2m 
M} - iM? 



4- • 2 7 4 

a J a, a,j — sin — a. 



2m — ajdj 



(20) 



9 u l 4- 

cos — aj 
2 J 



2m — djdj 



sin —\l 2m — d^dj dj 



cos Oj(m — d^dj) 



where N is the number of lattice's sites, e\ and ef. are 



Fermions' dispersions, 



= 2t(cos k x + cos k y + cos k z ) + \i (19) 
e% = — 2t(cosk x + cosky + cosk z ) + U — /i, 



where 9j = Q ■ r 3 and Q = (tt, 7T, it) is the antiferro- 
magnetic wave vector. In terms of the Bosc fields and 
keeping only the quadratic and quartic terms, the effec- 
tive Hamiltonian (EqfTT)) adopts the form 



l eff 



(21) 
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where 



h,2 = Jfn y ] ( a t a i + a j^ a j — a t a ^ ~ fflifflj) (22) 



Equation (|28|) shows that the Hartrcc-Fock parameter r 
renormalizes the exchange constant J. 

It is convenient to rewrite the Hamiltonian fEql28|) in 
momentum space representation: 



\ " („+ 



(ij) 



(«,' «,' "j <'j + afafafai 



(23) 



+ a^aidiCLj + aj~ajajCii ~ Adfdj'didj) 



and the terms without Bosc fields are dropped. 

The next step is to represent the Hamiltonian in the 
Hartree-Fock approximation. To this end one represents 
the product of two Bosc fields in the form 

a~laj = afa,j — < afcij > + < afa,j > (24) 

and neglects all terms (dfdj — < d[dj >) 2 in the four 
magnon interaction Hamiltonian. The result is 

a+a+a+aj ps - < a+a+ >< a+a 3 > 



afafa'j'ai 



— < ajdi >< ajaj > 
+ a+a+ < a+aj > +a+a; < a+d+ > 
afaididj ps — < afcti >< didj > 

+ a ,11 , < afcii > +afcti < a ,11 , > 



djdjdjdi 



2afa~j'a i aj 



~ — < a ■ Oj >< djdi > 

+ ajdi < a^aj > +a~j~aj < ajdi > 

PS — < afdj >< CLiCLj > 

- < a+a., >< a+dj > 

+ a+a+ < a q aj > +a, l a j < afaf > 

+ a^cii < djCij > +a~j~aj < afdi > 



(25) 



We assume that the matrix elements do not depend on 
the lattice's links and < dfdj >=< d%dj >. Then the 
Hartree-Fock approximation for the effective Hamilto- 
nian (Eq l21[) can be represented as a sum 



where 



h e ff ps tiHF — h c i + h q 



h cl = 6Jm 2 N(r- l) 2 , 



h q = Jmr (a + " 

(ij) 



d{ + a - dj — d i dj 



(26) 
(27) 

- didj) (28) 



and r is the Hartree-Fock parameter, to be determined 
sclf-consistcntly from the equation 



1 1 V- + 

= 1 > < dl dh > 



2m N 



(29) 



1 1 cos k x + cos k v + cos k z , , 

> < a7a j > 

2m TV 4^ 3 k - k 



h q = y~ K J [sa k d k - 7 fe (a+ a+ fc + a fc a_ fe )] (30) 



where 

e = 6Jmr (31) 
7fc = Jmr (cos fc^ + cos k y + cos k z ) 

To diagonalize the Hamiltonian one introduces new Bose 
field afe, by means of the transformation 

a fc = u fe a k + v k at k a~l = u k a k + v k a_ fc (32) 

where the coefficients of the transformation u k and v k 
are real functions of the wave vector k 



(33) 




The transformed Hamiltonian adopts the form 
h q =Y J {E k a+a k + El), 



with dispersion 



fees,. 



E k = 



and vacuum energy 



El 



47 



(34) 



(35) 



(36) 



For positive values of the Hartree-Fock parameter and all 
values of fc G 73, the dispersion is nonnegative E k > 0. It 
is equal to zero at k = (0,0,0) and k* = (±7r, ±7r, ±7r). 
Therefor, a^-Boson describes the two long-range exci- 
tations (magnons) in the spin system J35|. Near these 
vectors the dispersion adopts the form E k cx c s |k| and 
E k cx c s |k — k*| with spin- wave velocity c s = 

One can rewrite the equation for the Hartree-Fock pa- 
rameter (Eci l2T)|) in terms of the a k field 

r(T) - l + i-^lEv^Il + JMT)] 

k 

e k = cos k x + cos k y + cos k z (37) 
where n k is the Bose function of the a excitations 

1 



n k (T) 



e t — 1 



(38) 
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It is important to stress that the equation ([37)) can be 
obtained from the equation 

OF/dr = (39) 

where T is the free energy of a system with Hamiltonian 
h HF (EqUU) 

J- = 6Jm 2 (r-l) 2 + l^^ + ^ln (l - e"*) . 

k k 

(40) 

The sublattice magnetizations M A and M B for sub- 
lattice A (cos#; = 1) and sublattice B (cos#; = — 1) 
are defined by the equation (|20[) . It is evident that 
M A = —AI B , so that the total magnetization is zero. 
In terms of the Bose function of the a excitations 
they adopt the form 

M A (T) = -M B (T) = m - 1 Y, < > (41) 

k 

= m+ - - -— V , 3 [i +2n k (T)} 

At Neel temperature T/v the sublattice magnetization is 
zero M A (T N ) = -M B (T N ) = 0. From equation gTJ) 
and equation (|37[) rewritten at Neel temperature one ob- 
tains a system of equations which determines the Neel 
temperature 

k 

2m+1 = J-^ [l + 2n k (T N )} (42) 

To clarify the importance of the notion effective spin 
to one investigates the relationship between Neel tem- 
perature and to. The dependence of the dimensionlcss 
temperature Tjv / J on effective spin is depicted in fig- 
ure @. Decreasing the effective spin decreases the Neel 
temperature. The quantum critical value of the effective 
spin, the value at which Tn = 0, is m cr = 0.078. The 
effective spin decreases because the density of the doubly 
occupied states increases. The quantum critical point is 
a state with domination of the doubly occupied sites. 

Utilizing the dependence of the effective spin m and the 
exchange constant J/U on the parameter t/U (see figure 
[1]) one can obtain the dependence of the dimensionless 
temperature T^/U on the ration t/U. The phase dia- 
gram in plane of temperature Tjv/C/ and control param- 
eter t/U is depicted in figure @. The quantum critical 
value of the ratio is t/U = 0.9. 

To compare with experimental temperature-pressure 
curves one have to establish the relationship between 
hopping parameter t and pressure and between Coulomb 
repulsion U and pressure. The simplest assumption that 
U is a constant and t is a linear function of the pressure 
leads to a result which well reproduces the temperature- 
pressure phase diagram of CeRhIn§. But, to obtain the 
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FIG. 2: (Color online) The dependence of the dimensionless 
temperature Tn/J on the effective spin of the itinerant elec- 
tron m. The quantum critical value of the effective spin, the 
value at which Tat = 0, is m cr — 0.078 

experimental phase diagrams of CePd^Si^ or Celn^ one 
has to implement much more complicate fitting proce- 
dure. 




0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 



FIG. 3: (Color online) Phase diagram in plane of temperature 
Tn/U and control parameter t/U. The quantum critical value 
of the ratio is t/U = 0.9 



IV. PARAMAGNETIC PHASE 

When the system undergoes thermal (Tn > 0) or 
quantum (Tn = 0) transition to paramagnetic state, the 
magnon's dispersion opens a gap. This is a generic fea- 
ture of the second order phase transition. To describe it 
mathematically, one utilizes the modified spin-wave the- 
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ory proposed to describe 2D ferromagnetic [36J, |37[ and 
antiferromagnetic [HI, [39j| systems at finite temperature. 
The Takahashi's idea is to supplement the spin- wave the- 
ory with the constraint that the magnetization be zero. 
In the present paper we formulate, along the same line, 
a modified spin-wave theory of the paramagnetic phase. 

To enforce the magnetization on the two sublattices to 
be equal to zero in paramagnetic phase, one introduces 
parameter A, and the new Hamiltonian is obtained from 
the old one (Eq. [T7]) by adding a new term 



\^2(m-afa,i) 



(43) 



This modification leads to a modification of the Hamil- 
tonian fEq|3"U)). One obtains 



E 

k 



with 



6Jmr + A 



(44) 



(45) 



We implement the same calculations as above and arrive 
at a Hamiltonian which is modification of the Hamilto- 
nian (Eq[ 



keB 



with new dispersion 



(E k a + k 



(ik 



El 



Ei. = 



and new vacuum energy 



- 47 



4° 



(46) 



(47) 



(48) 



The free energy J 7 of a system with the modified Hamil- 
tonian reads 



T= 6Jm 2 (r- 1) 



i-V 

N 

k 



El 



-E 

k 



In 1 



(49) 

Then, one can obtain the system of equations for the 
Hartree-Fock parameter and the parameter A from the 
equations 



dT/dr = 0, dT/dX = 0. 



(50) 



The result is 
r(T) = 

2m + 1 = 



1 

4m 



(51) 



- 1 E 



3e 



12m N 

1 E — 



2Jmre\ 



e 



2h k {T)\ 



[l + 2n k (T)] 



where n k is the Bosc function of a excitation (Eq |38|) 
with new dispersion E k fEq|47l). 

It is convenient to represent the parameter A in the 
form 



A = QJmuK 



(52) 



introducing a new parameter k. Near to the zero 
wave vector the dispersion adopts the form E k oc 
v/cj |k| 2 + 36Jmr(2 k + k 2 ), where 36 Jmr (2k + k 2 ) is the 
gap of the magnon. It is zero below Neel temperature 
and increases when the temperature increases above Neel 
temperature. 

We implement the following procedure to calculate the 
Hartree-Fock parameter r and the parameter k. At tem- 
peratures below the Neel temperature k — and r(T/J) 
is obtained from equation ([3T| . At temperatures above 
the Neel temperature the functions r(T/J) and k(T/J) 
are solution of the system (fSTj) . The result is depicted in 
figures (gj) and ©. 

The figure ((4|) shows that the renormalization r at zero 
tempearture, due to the magnon-magnon interaction, is 
stronger when the system approaches quantum critical 
point (curve " a" ) , and it is insignificant for spin-localized 
systems (curve "d"). 



(a) t/U=0.9 




FIG. 4: (Color online)The dependence of the Hartree-Fock 
parameter on the dimensionless temperature T/J: (a) at 
quantum critical point t/U = 0.9 (m cr = 0.078); (b) at 
t/U = 0.425 (m = 0.175); (c) at t/U = 0.254 (m = 0.3); (d) 
at t/U = 0.09 (m = 0.5). The vertical dot lines correspond 
to the Neel temperatures Tn/J- 

The k parameter is a measure for the gap of the 
magnon in paramagnetic phase. It is zero at Neel tem- 
perature and increases when the temperature increases. 
The function k(T/ J) is depicted in figure ([S]) for different 
values of the ratio t/U. 

The figure ([5]) shows that the increase is faster when 
the system approaches quantum critical point (curve 
"a"), and it is weak for spin- localized systems (curve 
"d"). 
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FIG. 5: (Color online) The dependence of k parameter on 
the dimensionless temperature T/J: (a) at quantum critical 
point t/U = 0.9; (b) at t/U = 0.425; (c) at t/U = 0.254; (d) 
at t/U = 0.09. 



V. SPECIFIC HEAT 

Utilizing the above calculated functions r(T/J) and 
k(T/J) one can calculate the magnons' contribution to 
the specific heat of the system. By definition, the entropy 
is 



S 



dT 
dT 



dT dr 



dTdX 
~d\df 



dT 
dT 



(53) 



where T is the free energy of the system Eqs. (|40I49[) . Ow- 
ing to the equations (|39| and (|50| the first two terms in 
Eq. (1531) arc equal to zero and one obtains the customary 
formula for the entropy of a Bose system 



S 



-Y [(i 



+ n k ) ln(l + rife) - n k In rife] 



(54) 



where the dispersion E k Eq. (I35|) is used to define the Bose 
function Eq. ([38]) below the Neel temperature, and disper- 
sion iSfe Eq. f4"T)) above it. With entropy, as a function of 
temperature in mind, one can calculate the contribution 
of magnons to the specific heat. 



C V =T 



dS 
dT 



(55) 



The resultant curves C v (T/ J) , for different values of the 
parameter t/U, are depicted in figure ([6]). 

The function C V (T/J) has a maximum at Neel temper- 
ature. The figure © shows that the maximum is sup- 
pressed at quantum critical point (curve " a " ) ■ The max- 
imum is well observed exp er iment ally [32l. l40l |. and one 
can use it to determine the Neel temperature. The ex- 
perimental measurements of the specific heat of CeRhln^ 
for different pressures [32| show that at TV, C v /T has a 
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0,8 1,0 1,2 1,4 
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FIG. 6: (Color online) Contribution of magnons' fluctuations 
to the specific heat: (a) at quantum critical point t/U = 
0.9 (m cr = 0.078; (b) at t/U = 0.425 (m = 0.175); (c) at 
t/U = 0.254 (m = 0.3); (d) at t/U = 0.09 (m = 0.5). 



very sharp peak at ambient pressure. With increasing 
pressure the magnetic anomaly remains well defined but 
maximum decreases. This phenomena is well described 
theoretically in the present paper Fig.©. The existing 
correspondence between specific heat and derivative of 
the resistivity suggests a confidence that magnon's fluc- 
tuations are important and for the transport properties 
of the itinerant antifcrromagnets . 



VI. SUMMARY 

In this paper an itinerant antifcrromagnets were stud- 
ied. Varying the ratio t/U, where t is the hopping pa- 
rameter and U is the Coulomb repulsion, the system 
was investigated between a state with localized spins 
(t/U < 0.09) and quantum critical point t/U = 0.9 at 
which Neel temperature is zero. The evolution of the 
magnons' fluctuations in antifcrromagnctic and param- 
agnetic phases was studied by means of the renormalized 
spin-wave theory and modified spin-wave theory. The 
renormalized spin-wave theory includes a parameter r 
which accounts for magnon-magnon interaction. The sys- 
tem of equations for the Neel temperature (|4"2j) rewritten 
for 2D system has the only solution TV = 0. This shows 
that present method of calculation is in accordance with 
Mcrmin- Wagner theorem, which is our theoretical crite- 
ria for adequate account of magnons's fluctuations. The 
modified spin-wave theory involves in a natural way the 
gap of the magnon in paramagnetic phase. 

The effective spin is the most important factor which 
determines the quantum criticality. To figure out this 
one has to map the low-energy excitations of the half- 
filled Hubbard model onto an effective spin- 1/2 Heisen- 
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berg model [24|. The dimcnsionless Neel temperature 
T N /J for this model is T N / J = 1.387. The result shows 
that the Neel temperature of the spin- 1/2 antifcrromag- 
nets is nonzero for all values of the exch ang e constant 
including weak coupling U/t < < 1 regime [24| . This is in 
contrast to the result in the preset paper where the an- 
tifcrromagnctism is destroyed at U/t < l.ll(t/U > 0.9). 
The mapping of the Hubbard model onto an effective 
Heisenberg model with effective spin- to < 1/2 is a crucial 
step towards the understanding of quantum criticality. 

The critical value of the effective spin to = 0.078 only 
depends on the effective Heisenberg model. The effec- 
tive exchange constant J and the effective spin to are 
calculated in one fermion-loop approximation. It is nei- 
ther strong coupling approximation nor weak coupling 
approximation. At half filling the chemical potential is 
fixed fi = U/2 and the dispersions' Eq. (TT5]) dependence 
on the parameter t/U is nontrivial. 

To compare with experimental results the phase dia- 
gram in plane of T/U and control parameter t/U was ob- 
tained. To compare with results of the numerical calcu- 
lations one has to convert the diagram Fig.Q into phase 
diagram in plane of T/t and U/t. The result is shown in 
figure 0. 



0,8 
0,7 



^^^present paper 
-•- DMFT 

O DCA 
— A — DrA 




FIG. 7: (color online) Phase diagram in plane of temperature 
Tn /t and control parameter U/t. The solid black line is the 
phase diagram obtained in the present paper. The red circles, 
the blue open diamonds and the magenta triangles are the 
results in (DMFT), (DCA) and (DrA) respectively. 



The solid black line is the phase diagram obtained in the 
present paper. The red circles, the blue open diamonds 
and the magenta triangles are the results obtaine d by 
means of the dynamical mean field theo ry ( DMFT) [28| . 
dynamical cluster approximation (DCA) [28[ and dynam- 
ical vertex approximation (DrA) (30j , respectively. 



One of the most important point in the Hubbard model 
is the maximum of the Neel temperature Tjv /t as a func- 
tion of U/t. It indicates the crossover from itinerant mag- 
netism to the magnetism of the localized spins. The Neel 
temperature of itinerant systems increases, when U /t in- 
creases, as a result of the increasing of the effective spin. 
The Neel temperature of the spin- 1/2 Heisenberg model 
of localized electrons is = 1.387J with exchange con- 
stant which decreases when U/t increases. In the present 
paper the maximum is at U/t = 10.6, in the DMFT 
at U/t = 10 and in DrA at U/t = 9.8. The overall 
agreement of the results for the Neel temperature in the 
present paper with the results obtained in DMFT are sat- 
isfactory for all values of U/t except for lowest one. The 
non-local corrections, accounted for in DrA, reduce the 
Neel temperature Tjy/t versus temperatures in DMFT 
and present paper in whole phase diagram. 

The exchange constant J of the effective model Eq. 
(|17p is calculated in the limit when the frequency and the 
wave vector are small. The calculations can be improved 
employing an effective Heisenberg theory with exchange 
constant J(k) which depends on the wave vector k. 

Another important point in the Hubbard model is the 
quantum critical point (QCP) Tn = 0. The phase di- 
agram Fig.Q shows that it is reached at U/t = 1.11. 
The figure also shows that all numerical calculations are 
implemented at control parameters U/t > 3.674, which 
is far from the QCP. This explains the absence of com- 
ments on the quantum criticality. On the other hand an 
investigation of the Hubbard model without explicit affir- 
mation of existing or nonexisting of QCP in the Hubbard 
model is not compleat. One can extrapolate the curves 
in DMFT, DCA and DrA down to zero temperature. A 
non-zero QCP emerges in these approaches. 

Alternatively, the Hubbard model is studied by a map- 
ping on a spin- 1/2 generalized Heisenberg model with 
higher order terms in a form of long-range or ring ex- 
change [4l|. An effective Heisenberg model with spin- 
1/2 is considered which can not describe the quantum 
criticality despite the fact that parameter t/U is close to 
the quantum critical point. This is because the quantum 
critical behavior depends decisively on the effective spin 
of the itinerant electron which is smaller then 1/2 near 
the quantum critical point (see figure ©). 

Neel order and quantum critical point in the Hubbard 
model arc also studied by means of a spin-charge rotating 
reference frame approach [42j , which explicitly factorizes 
the charge and spin contribution to the electron opera- 
tor. The effective constants are calculated by means of 
Hartree-Fock approximation of the fermion interaction. 
This explains the over estimation of the critical value of 
the control parameter U/t = 0.676(t/U = 1.479) com- 
pare with the result in the present paper t/U = 0.9, 
where the Coulomb repulsion is treated exactly. 
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